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of the Fermilab Tevatron (^/s = 2.0 TeV) is revisited. For Higgs masses below 
140 GeV we demonstrate that the production cross section times branching ratio 
for H — > bb decays yields a significant number of events and that this mode 
is competitive with and complementary to the searches using pp — > WH, ZH 
associated production. For higher mass Higgs bosons the H — > W + W~ decays 
are more difficult but have the potential to provide a few spectacular events. 

One of the untested features of the elegant and highly successful Standard Model (SM) of 
electroweak interactions is the Higgs mechanism for spontaneous electroweak symmetry breaking. 
The Higgs sector of the SM consists of one doublet of complex scalar fields which is used to break 
the SU(2)xU(l) electro-weak gauge symmetry to U(l) electromagnetism. This mechanism gives 
rise to the masses of the W and Z bosons and is also used to give masses to the charged leptons and 



quarks. An additional physical neutral scalar particle is also produced, the undiscovered Higgs 
boson. All properties of this boson are predicted by the SM except for its mass. Most extensions 
to the SM's electroweak symmetry breaking mechanism predict more than one physical scalar 
state, however one of these states tends to have properties very similar to the SM Higgs boson. 

Direct collider searches have so far ruled out a SM Higgs boson with a mass less than « 
107 GeV DU at the 95% confidence level. The final run of the Large Electron Positron collider 
(LEP) at CERN this year should be able to extend the exclusion limit or provide a 5<r discovery 
up to about 114 GeV. The next CERN collider, after LEP, is the 14 TeV proton-proton Large 
Hadron Collider (LHC), which will have the capability to discover a SM Higgs boson of nearly 
any massQ in multiple search channels. However, the LHC is not scheduled to begin taking data 
until 2005 at the earliest. In the meantime, the upcoming Run II of Fermilab's 2 TeV proton- 
antiproton collider, the Tevatron, will have an opportunity to discover a Higgs boson above the 
LEP exclusion limit and below 180 GeV. In this Letter we will concentrate on the SM case over 
this mass range. 

The SM Higgs boson couples at lowest order to all SM particles with a strength proportional 
to their mass. For the massless photon and gluon a coupling is induced at higher order via loops of 
heavy particles. Therefore there are many ways for the Higgs boson to decay and also many ways 
to produce it at a hadron collider. For the mass range of interest here, 100 ^ Mh ~ 200 GeV, 
the Higgs boson decays predominantly to either a pair of 6-quarks (bb) for a mass less than 135 
GeV, or to four fermions via real or virtual W gauge bosons for a mass above 135 GeV. These 
decays are referred to as H — > bb and H — > W + W~ respectively, and at 135 GeV they have 
approximately equal branching ratios (43% to 40% respectively). At this mass, the remainder of 
the decay branching ratio is primarily into r lepton, c quark and gluon pairs. 

Cross sections for various Higgs boson production modes at the Tevatron are shown in 
Fig. [l|(a) [0]. The dominant mode is via gluon-gluon fusion, gg — > H, which proceeds via a 



In the SM Mjj > 1 TeV is ruled out by unitarity constraints. 
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top quark loop. The next largest mode is Higgs boson associated production with a W or Z 
boson. At a still lower cross-section is Higgs boson production via bremsstrahlung off a top 
quark pair. The current searches planned for the Tevatron Run II concentrate on the domi- 
nant production modes: WH,ZH associated production with decays H — ► bb ||, and gluon 
fusion production gg — > H with subsequent decay H — > W + W~ [|Tj. While detailed studies with 
detector simulation and neural net analysis show promise, they do indicate the need for large 
integrated luminosity, and even with 30 fb _1 these searches have difficulty in the mass regions 
130-150 GeV and above 180 GeV. In this letter we investigate the feasibility of a search in tiH 
associated production, first proposed in Refs. || in the context of the Superconducting Super- 
collider, and later examined in more detail for the LHC in particular; see || and 0. Decays to 
both H — > bb and H — > W + W~ are considered, with a varying number of tagged b quarks and 
additional jets in the final state and at least one charged leptonP] to help discriminate against 
the large QCD backgrounds. We believe that, in spite of the low cross section, a search for the 
SM Higgs boson in association with a top quark pair in the Tevatron Run II data set (15 fb _1 
integrated luminosity) is likely to be comparable to the early search and evidence for the top 
quark results from CDF in Run I 

We have calculated the signal at the parton level for pp collisions at ^/s = 2.0 TeV using exact 
tree-level matrix elements generated by madgraph |§ and comphep JkJ, and NLO corrected 



decay rates of the Higgs boson via hdecay ||11|| . We take the K-factor to be 1.33: the ratio of the 
NLO ti cross section to the leading order value. Our estimate agrees well with a calculation in 



the effective Higgs approximation [12], which estimates the K-factor to be 1.2 to 1.5 depending 



2 Top quarks decay predominantly to a IT-boson and a 6-quark. The IT-boson subsequently decays 2/3 
of the time to a pair of quarks (jets) and 1/9 of the time into each of ev e , \xv^ and tv t . The neutrinos 
are not observed and the r-lepton further decays 2/3 of the time to a jet plus neutrinos and 1/3 of the 
time to e or /x and more neutrinos. Thus, about 26% of the time a TT-boson decays into a taggable 
charged lepton (e or [i). 
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FIG. 1. For pp collisions at y/s = 2.0 TeV: (a) cross section (pb) for various Higgs boson production 
modes (b) the tiH production cross section times branching ratio and K- factor (fb) for H — > bb 
(solid) and H -» W + W~ (dashed). 



on the Higgs boson mass and other uncertainties. CTEQ4L |T3j parton distribution functions are 
used, and both the factorization and renormalization scales are taken as the top quark mass, m t . 
The cross section times Higgs branching ratio times K-factor used throughout the rest of this 
letter are shown in Fig. 0(b). For H — > bb the cross section times branching ratio falls steeply 
and becomes smaller than 1 fb for Higgs masses above 140 GeV. In contrast, the cross section 
times branching ratio for H — > W + W~ is a broad bump above 1 fb between 125 and 190 GeV. 
Thus, a search strategy developed for a 160 GeV Higgs would likely work over a much wider 
mass range, covering some of the unobservable regions of the current searches. 

For lower Higgs boson masses (< 140 GeV), the most interesting decay mode of the Higgs 
is H — > bb and the final state of the hard interaction is W + W~bbbb. For higher mass Higgs 
bosons, H — > W + W~ is the dominant decay mode and the final state signature is instead 
W + W~W + W~bb. For both modes, each W-boson will decay either hadronically into a quark 
and anti-quark or leptonically into a charged lepton and neutrino. Backgrounds to these modes 
consist of resonant production, e.g. irreducible tiZ events with Z — > bb or Z — > £ + £~; continuum 
tibb or ticc production; and reducible backgrounds such as ti + jets, where the additional jets 
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TABLE I. Inclusive cross sections for the major physics and reducible backgrounds for pp collisions 
at y/s = 2.0 TeV. Continuum production of final states involving heavy quarks or weak bosons plus 
additional jets include a px > 20 GeV cut on the additional jets, but no cut on the decays of heavy 
states. Typical uncertainties on these estimates are on the order of ±50%. 
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are light quarks or gluons but may be misidentified as b quarks. Table | itemizes the major 
backgrounds and gives their cross sections. All other SM backgrounds are negligible. Clearly, 
the largest background arises from ti + jets. Thus, the search for tiH will critically depend on 
the detector ability to tag b quark jets and simultaneously to suppress mistagging of non-6 jets. 

We have calculated the parton level cross section for ti + jj events using the exact tree-level 
matrix elements 0,0], which contain both collinear and soft singularities for the additional jets. 
The collinear singularities are removed by imposing the requirement that the final-state partons 
be well-separated in space, AR(jj) > 0.4. This corresponds to the experimental requirement that 
the observed jets be similarly separated. f\ The severity of the soft singularity for a given pt cut 
may be estimated by examining the ratio of the cross sections for ti+jets production to that for ti, 
the latter of which may be regarded as an inclusive rate; the overall rate can later be normalized 
to NLO calculations. For the Tevatron we find Ottjjl&tt ~ 1/7, including a requirement of 
Pt > 20 GeV for the additional jets, indicating that the calculation is largely perturbative. This 
is in line with general expectations that each additional gluon jet of pt ~ 20 GeV would multiply 



5 Ai? is a separation in the space of detector coordinates azimuth and pseudorapidity, AR(jj) 



y/A(j)' 2 + Al]' 2 . 
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the inclusive rate by a factor « 1/2 — 1/3. One may also use an exponentiation approximation 
for the soft gluons [[OJ, which yields a result 30% lower than the value in Table |. 

We have compared the results of our matrix element calculation with Pythia 6.115 |TB| , and 
while there are differences in cross section as a function of jet p™ in , the differences at Tevatron 
energies are largely within the uncertainties of the matrix element calculations. Ultimately, any 
uncertainty in the rate for additional jets in ti events is not a concern as the ti + jets sample in 
Run II will be large enough to calibrate the Monte Carlos. Furthermore, by retaining the larger 
estimate we make this analysis more conservative. 

These uncertainties aside, it is clear that a major background could come from the mistagging 
of non-6 jets. Both the efficiency for b jet tagging and non-6 jet suppression will thus be extremely 
important in this search. Our studies of the upgraded Tevatron collider detectors using Run I 
algorithms predict an average b jet tagging efficiency of about 60% and a c quark mistag rate 
of about 25%, while at the same time providing at least a factor of 500 suppression of light 
quark and gluon jet mistags. We estimate that implementation of 3-D vertexing algorithms will 
improve the b tagging efficiency to at ~ 70% while leaving the mistagging rates unchanged. We 
use the 3-D values in our analysis (70%, 25% and 0.2%). 

For a Higgs boson of mass Mu < 140 GeV, where the W + W~bbbb signature is preferred, a 
70% b tagging efficiency translates to a 92% probability to tag > 2 b jets, 65% for > 3 tags, 
and 24% to tag all four. As an example, consider the case of Mh = 120 GeV, with one W 
decaying leptonically and one decaying hadronically, which has much better known backgrounds 
than the larger all-hadronic W decay channel. The cross section times branching ratio is 2.0 fb. 
For three b tags, 15 fb -1 of data would contain ~ 19 signal events and about 140 background 
events. Requiring four b tags would yield a signal sample of approximately 7 events, on top of a 
background of about 40 events. These counts are before any top reconstruction or Higgs boson 
mass binning. 

To extract these few signal events we study their reconstruction by performing a generator 
level Monte Carlo simulation using Pythia 6.115, which is in good agreement with our parton- 
level calculations for tibb production using a jet separation corresponding to a parton-level cone 
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cut of AR > 0.7. Event selection is determined using a parameterized detector simulation. 
For reasonable rejection of the ti + jets background, we require at least one isolated lepton 
and 4 jets with pr > 15 GeV, at least 3 of which are b tagged, and two additional jets with 
Pt > 10 GeV as well as a missing Et > 15 GeV. The tibb (ttcc) background differs from signal 
primarily in that the invariant mass of the b (c) quarks from gluon splitting will be quite low, 
whereas the b quarks from the Higgs boson will have a much higher invariant mass. Prior to top 
reconstruction, however, in any given signal (background) event we do not know which of the 
four b jet candidates are from the Higgs boson (gluon). Thus, for each event we form invariant 
masses of all combinations of these b jet candidates and order them. We find there is considerable 
separation between signal and background in the fourth-highest of these ordered dijet masses, 
and cut at 60 GeV in this distribution, see Fig. @(a). Incorrect selection of a jet from the hadronic 
W decay is taken into account. 

In Run I, ti reconstruction efficiencies of 60% were achieved by CDF for double 6-tagged 



events [fT? ]. Studies indicate that this can be enhanced by improvements in jet-parton assignments 
resulting from better tracking and energy corrections expected in Run II. Further improvements 
may even be possible via detector enhancements to b flavor tagging. However, the efficiency for 
tiH may be lower than for ti, but we do not make the distinction here. In this paper we will 
use an optimistic 70% efficiency but note that results are only slightly worse for 60%. Although 
this is somewhat higher than that of the LHC studies, 0, the jet multiplicity at the Tevatron is 
significantly lower than at the LHC, and as the tiH events would be produced closer to threshold 
at the Tevatron, they would be better measured than at the LHC. Using 70% efficiency we plot in 
Fig. H](b) the invariant mass distribution of the additional jets after top quark pair reconstruction. 
Here the signal exhibits a significant peak above the background. For Mh = 120 GeV, this 
would correspond to a 2.8a observation for one experiment using Poisson statistics converted 
to a Gaussian equivalent, or 4.1a for two experiments. Increasing the integrated luminosity to 
20 fb _1 increases the significances to 3.3a /4.7a, respectively. 

For higher Higgs boson masses, M# > 140 GeV, the dominant signature is W + W~W + W~bb. 
One possible strategy is to demand > 3 charged leptons, which is one of the smaller components 
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FIG. 2. (a) The tiH signal (solid) and ti + jets (including tibb, ticc) background (dotted) distri- 
butions for the fourth-highest mass pair of 6-tag candidates. This includes cases where a jet from the 
hadronic W decay is misidentified as a b jet. We cut on this distribution at rribb > 60 GeV. (b) Invariant 
mass distribution of 6-tag candidates after top quark pair reconstruction with efficiency e rec = 70%. For 
M H = 120 GeV and 15 fb" 1 this corresponds to a 2.8a significance for one experiment, or 4.1cr for two. 

of this decay mode at ~ 6% of the total rate. For M# = 160 GeV and 15 fb _1 of data for one 
experiment, this translates to only 3.0 events of the 45 total tiH events produced. Adding a 
requirement to see > 1 b tag and minimal cuts, the signal sample is only 2 events. The expected 
background, however, is much less than one event, so there is the potential to observe one or two 
spectacular events per experiment^. This estimate holds approximately over the Higgs boson 

4 Any background from ttW may be eliminated by requiring the observation of extra hadronic activity 
from the fourth signal FF-boson. 
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mass range 140-170 GeV. 

The signal sample containing exactly two charged leptons is larger, by a factor of 3.6, but 
there is much more background. Requiring the two leptons to have the same charge reduces 
the background dramatically, to of order 2 events if additional jets are required to be observed. 
Further mass constraints may make this much smaller than one. However, since only 1/3 of the 
two lepton signal sample will have same-sign leptons, one expects another 3-4 spectacular events 
per experiment. 

Despite the low signal rate for ttH associated production events, there is considerable po- 
tential in this channel to observe a SM Higgs boson, or to set exclusion limits on its mass. We 
have not examined cases of non-SM Higgs bosons, but the outlook is optimistic as production 
cross sections times branching ratio are often enhanced relative to the SM, e.g. in SUSY or 
topcolor models. Our analysis suggests an optimistic search strategy for Higgs bosons up to 
about 140 GeV in the H — > bb channel, and some promise in the H — > W + W~ channel for 
the mass range 140-170 GeV. A thorough understanding of and improvements to the top quark 
reconstruction efficiency and the identification of b, c and light quark jets are very important 
for this search. Increasing the integrated luminosity of the Tevatron by a factor of three would 
increase the significance of this search to 5a per experiment. 
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